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ABSTRACT 

This paper discusses the evolution of the correlation between galaxy star formation rates 
(SFRs) and stellar mass (A/,) over the last ^10 Gyrs, particularly focusing on its environ- 
mental dependence. We first present the mid-infrared (MIR) properties of the Ha-selected 
galaxies in a rich cluster C10939H-4713 at z = 0.4. We use wide-field Spitzer/MIPS24/im 
data to show that the optically red Ha emitters, which are most prevalent in group-scale envi- 
ronments, tend to have higher SFRs and stronger dust extinction than the majority population 
of blue Ha sources. With a MIR stacking analysis, we also find that the median SFR of Ha 
emitters increases in higher-density environment at z = 0.4, and this trend is confirmed for 
both red and blue galaxies. The trend becomes much less significant when we compare their 
specific SFR (SSFR), although we find that there still remains a weak, but positive correlation 
between SSFR and galaxy number density. We then discuss the environmental dependence 
of the SFR versus A/, relation for star-forming galaxies since z ^ 2, by compiling com- 
parable Ha-selected galaxy samples in both distant cluster environments (from MAHALO- 
Subaru) and field environments (from HiZELS). Based on these large, narrow-band selected 
Ha emitter samples, we find that the SFR-Af* relation does not show detectable environmen- 
tal dependence since z ^ 2. Our data also suggest that the SSFR of Ha-selected galaxies (at 
log(A/»/A/0) =10) evolves as (1 + z)^, and this evolutionary speed seems to be independent 
of the environment. Therefore we suggest that the primary star-formation quenching mecha- 
nism in cluster environments is a fast-acting process at any time in the history of the Universe 
since z ~ 2. 

Key words: galaxies: clusters: individual: CI 0939H-47 13 — galaxies: evolution — large-scale 
structure of Universe. 



1 INTRODUCTION 

Galaxy formation and evolution is strongly dependent on envi- 
ronment. In the local Universe, galaxies in cluste r environments 
are m ostly passive (red), early-type galaxies (e.g. lDresslei|[T980l : 
lOoto et al. 2003), and there is a clear trend that the star forma- 
tion activity of galaxies tends to be lower in high-density environ- 
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ment than low-density fields (e.g. Lewis et al]|2002l : iGomez et al.l 
l2003l : iTanaka et alj|2004l : iBalogh et al.ll2004h . Therefore it is be- 
lieved that the star formation activity of galaxies is affected by 
their surrounding environment during the course of the cluster or 
group assembly process. The environmental trends are also seen in 
the distant Univers e (e.g. lKodama et al .l200ll : |Postman et alj200l 
lOuadri et al]|2012h . and some recent studies suggest that the star 
formation-density relation may be reverse d at z > 1 jElbaz et al.l 
I2OO7. : .Cooper et al...200a: .Ideue et al..,2009t iTran et aLll2010l) . This 
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"reversal" of the star for mation-den sity relation in the early Uni- 
verse is still under debate JPatel et ai., 2009: 2011), because the re- 
sults could be highly unce rtain depending on the sample definitions 
or selection method (e.g. iFeruglio et"alll2010l ISobral et ailboill : 
IPopesso et al.ll20l"lh . Nevertheless, some recent studies which fo- 
cus on individual galaxy clusters indeed find a hint that a fraction of 
galaxies in distant {z ~ 1) cluster outskirts or inter mediate-density 
enviroimients are showing boost ed activity (e.g. Marcillac et al 



2007| ; |Poggianti et alj200a 2009: lKovama et alj201Cl . lGeach et al 



20111) . suggesting an accelerated galaxy evolution at the site of ac- 
tive cluster assembly. 

A growing number of studies have revealed important roles 
of the group-scale environment or in-falling regions around 
rich clust ers, by s tudying e.g. galaxy colours or morph olo- 
gies (e.g. [ Kodama et al. 2001; Wilman et al. 2008; Balogh et~aLl 
I2OI ih . In addition, optical emission-line surveys or MIR-FIR 
observations of distant clusters have brought some new in- 
sights into the (obscured) nature of star forming galaxies in 
cluster environment. A prominent example is dusty red galax- 
ies, which are reported to populate the outsk irts of rich galaxy 
clusters out to z ~ 1 (e.g. Koyama e t al. 200^ 2010 ; 
see also Iwolf. Gray, & MeisenheimeJ 2005 : [Geach et atl I2OO6I; 



Verdugo, Ziegler. & GerkenI I2OO8I; IWolfetal.l |2009| ; iTranetall 
20091) . These studies suggest that such dusty red galaxies are a key 



population for understanding the physics of environmental effects. 

In this respect, wide-field HaA6563 emission-line surveys of 
distant galaxy clusters are a powerful method to pinpoint the loca- 
tion of this key population. This is not only because the Ha line 
is less affected by dust extinction compared to star-formation in- 
dicators at rest-frame ultra-violet wavelengths (e.g. [Oll]A3727), 
but also because emission line surveys with narrow-band (NB) fil- 
ters allow us to effectively pick out galaxies from a narrow red- 
shift slice. The latter is particularly important for cluster studi es be- 
cause the effect of contamination could be a concern. iKoyama et al.l 
( I2OI l[) performed a wide-field Ha emission-line survey of a rich 
cluster, C10939-I-4713 {z = 0.41), and find a strong concentra- 
tion of optically red star-forming galaxies in the group- scale en- 
vironm ent around the cluster. While we argued in Koya ma et al.l 
1 I20T l|) that the excess of the red star-forming galaxies suggests an 
enhancement of dust-obscured star formation in the group environ- 
ment, a firm conclusion still await a direct measurement of dust- 
enshrouded star-formation in those galaxies, because even Ha lines 
are repo rted to be heavily extinguished in extreme l y dusty galax- 
ies (e .g. IPoggianti & Wul I2OO0I ; lOeach et alJ[20o3 : iKovama eTai] 
I2OIOI) . Therefore, the first goal of this paper is to directly unveil 
the nature of this red star-forming population in distant group envi- 
ronments using mid-infrared (MIR) observations. 

Another important parameter that drives galaxy evolution is 
the stellar mass (M«) of galaxies. A correlation has been claimed 
between galaxy SFR an d M, for st ar-forming galaxies in the lo- 
cal Universe (e.g. Brinchmann et al. 2004; Peng et al. 201 0), as 
well as in the distant Universe out to z > 2 (e.g. |Noeske et al.l 
20071; iDaddi et alj|2007l: ISantini et aljlioogi : iKaiisawa et alj|2010l ; 



Bauer et al.l 1201 ll ; IWhitaker et al.l l20I2h . A growing number of 



studies are now investigating various aspects of this relation; 
e.g. the origin of its scatter or its morphological dependence 



jWuyts et aUHoTTI ; ISalmi et alJl2oT^ . An interesting implication 
from some detailed studies of the local SFR-Af, relation is the "in- 
dependence" of this relation with environment; IPeng et al.l JioToh 
studied local star-forming galaxies drawn from SDSS to show an 
excellent agreement in the SFR-A/, sequence between low-density 
and high-density environments. They argue that the environment 



does change the star-forming galaxy fraction, but that it has very 
little impact on th e SFR-Af, relation of t hose galaxies that are star- 
forming (see also lWiiesinghe et al.ll2012h . 

An observational challenge is to test this universality of SFR- 
M« relation in the distant Universe. Any local relation may not 
necessarily be applicable for distant galaxies at z > 1, where 
the average star- formation activity is abou t an orde r of magni- 
tude higher (e.g. 



l2013l ; IStottetaLl 



Hopkins & Beacomlbood see also ISobral et al.l 
20131) . Unfortunately, constructing a large, uni- 



formly selected galaxy sample at such high redshifts is still chal- 
lenging, which prohibits us from understanding the environmen- 
tal impacts on the SFR~A/, relation in the distant Universe. Some 
earlier works have attempted to identify the environmental depen- 
dence (or its absence) of the SFR-A/, relation out to z = 1 



t Vulcaniet al. 2010 


. Li et al. 201 1 ; McGee et al. 201 iLlTvler et al.1 


2011; Muzzin et al. 


2012), or out to z = 2 dTanaka et al. 


201C|; 


Tanaka et al.^ .201 1.; 


Griitzbauch et al.ll201 ll: iKoyama et al. 


2OI3I). 



but a full consensus on the envirormiental impacts on the SFR ver- 
sus A4^, relation has not yet been reached, because of the different 
sample selection and/or different environment definitions. There- 
fore, the second goal of this paper is to test the environmental 
dependence of the SFR-Af, relation out to 2 ~ 2 for the first 
time based on the purely Ha-selected star-forming galaxy sam- 
ples established in our recent two narrow-band Ha survey projects; 
MAHALO-Subaru for clusters and HiZELS for general fields (see 
details in § 2.2). 

The structure of this paper is as follows. In § 2, we present 
the wide-field Subaru and Spitzer/MlPS 24fim data of CI 0939, as 
well as the Ha emitter samples at different redshifts and in dif- 
ferent environments from our previous studies. We also describe 
how we estimate fundamental physical quantities such as stellar 
masses and star formation rates of galaxies. The main results and 
discussion of the paper are described in § 3. We first discuss MIR 
properties of Ha emitters around CI 0939, particularly focusing on 
the nature of the red Ha sources (§ 3.1-§ 3.2), and then we dis- 
cuss the environmental dependence of the SFR versus Af* relation 
since 2 ~ 2 using all the Ha selected galaxies (§ 3.3-§ 3.6). Fi- 
nally, our conclusions are given in § 4. Throughout the paper, we 
adopt the standard cosmology with JIm = 0.3, Qa = .7, and 
Ho = 70kms"^Mpc~\ and we assume a ISalpeted ( IT95I initial 
mass function (IMF). Magnitudes are all given in the AB system. 



2 DATA AND ANALYSIS 
2.1 The 01 0939 cluster 

2.1.1 Subaru data 

The CI 0939 cluster at z = .41 is one of the best-studied clusters at 
intermediate redshifts (e.g. Dressier c& Gunn 1992; Dressier et alj 
I994; Stanford, Eisenhardt, & Dickinson 1995; Small et al. 199|; 
Sato & Martini 120061 : iDressler et al.l 120091 lOemler et al.i |2009|) . 
In addition to these studies which focus on cluster central 
wide-field (~30'x30') optical broad-band and 



regions, very 

narrow-band i maging surveys of this clu ster have been made by 



^rnagu _ 

iKodarna et al f2001i) and . Koyama et alj [2011 ) using Suprime- 
Cam jMiyazaki et al.l I2OO2I) on the Subaru Telescope jive et al.1 
l2004h . The details of the se Subaru data have already been presente d 
in our previous papers dKodama et al.ll200ll : iKovama et al]|201lh . 
and so here we provide only a quick su mmary of th e data. The 
broad-band data are analysed by Kodama et alj j200ll) . They dis- 
covered a 10 Mpc-scale filamentary large-scale structure around the 
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cluster based on the photometric redshift (photo-z) techn ique. The 
Hq emitter search of this field is made bv .Kovama et al. ( 201 ih . us- 
ing the narrow- band filter NB921 (Ac=9180A) on Suprime-Cam, 
who identified >400 Hq emitting galaxies around the cluster. In 
this paper, we use both the Hq emitter and the photo- 2: selected 
cluster member catalogues presented in these studies. The photom- 
etry of the sources is pe rformed with the SEXTRACTOR software 
teertin & Amoutslll996h . We primarily use MAG.AUTO as the to- 
tal magnitudes for measuring physical quantities such as Ad, or 
SFR, while we use 3" aperture (corresponding to 16 kpc) photom- 
etry for measuring galaxy colours. 



2. 1.2 Spitzer MIPS 24fiin data 

We retrieve the wide-field MIPS (Rieke et al. 2004) 24^m scan 
data of the CI 0939 field from the Spitzer Science Archive. The data 
covers a large part of our Suprime-Cam field of view (from cluster 
core to surrounding groups), so are well-suited for studying dust- 
obscured star formation activity around the cluster. The data was re- 
duced from the Basic Calibrated Data stage (provided by the Spitzer 
Science Centre) using the MOPEX software following the proce- 
dure outlined in Geach et al. (2006). Source ex traction was per- 
formed using SExtractor dBertin & Amoutslll996b , with the criteria 
that a source consist of at least three contiguous pixels (each pixel 
is 2.5" square) at >2a above the background. We measure 16" di- 
ameter aperture fluxes, corresponding to ~3 times the FWHM of 
the point-spread function (PSF; FWHM of 5" at 24/im). Using a 
curve-of-growth analysis on bright isolated point sources, and our 
completeness simulations, we find that 16" apertures recover ~75 
per cent of the total flux, and we therefore correct the resulting 
fluxes by a factor of 1.33 to yield the total 24^m fluxes. 

The 24^m catalogue contains 886 sources down to 200/iJy 
(~5cr limit) within the Suprime-Cam field. We perform cross- 
matching between the 24^m sources and our optical sources to con- 
struct a MlPS-detected cluster member catalogue. We search for 
optical sources within 3" radius from each 24^m source, and we 
select the nearest source as its optical counterpart. As a result, we 
find 162 MIPS sources which are likely to be associated with clus- 
ter member galaxies (photo-z members or Hq emitters), although 
due to the poor spatial resolution of the MIPS24/im data, it is some- 
times difficult to identify the counterparts correctly. In order to 
avoid such ambiguous detections, we check all the sources by eye 
and exclude some heavily blended sources. Only 10 sources (6%) 
are excluded in this process (and they seem to be just normal galax- 
ies randomly distributed on sky), so the exclusion of these sources 
does not affect our results at all. Overall, the MlPS-detected mem- 
ber catalogue contains total 152 sources, among which 33 sources 
are Hq emitters. The relatively small number of the Hq detected 
sources is not surprising, because the photo-2 selected members 
include galaxies over a relatively broad range in reds hift (we apply 
030<2;phot<0.45 for photo-z member criteria; see iKovama et al.l 
l20Ilh . while the Hq emitters are considered to be secure cluster 
members located within the narrow redshift slice at the cluster's 
redshift (0.39<z<0.41). 



2.1.3 MIPS stacking analysis 

We note that the limiting flux of the 24^m data corresponds 
to SFR~8A/0 yr^^, while our Hq survey reaches down to 
SFR<Cli\/0 yr^^ (without dust extinction correction). Indeed, the 
number of Hq emitters which are individually detected at 24pm is 



not large, because of the limited depths of the 24^m data. We there- 
fore apply a 24/im stacking analysis for the Hq emitters around 
CI 0939 to study the general properties of the fainter sources, by 
dividing the sample into subsamples selected by colour or environ- 
ment (see § 3.2, § 3.3). We exploit median stacking at the positions 
of Hq emitters down to an dust-uncorrected SFR of O.25i\f0 yr^^ 
(which corresponds to a 5a detection in Hq). Although our conclu- 
sion does not change even if we apply an average stacking, median 
stacking is preferred because it can minimize the effects from ex- 
ceptionally luminous sources, as well as those from some luminous 
nearby sources. 

The photometry on the MIPS stacked image is performed in 
the same way as for individual sources; i.e. 16" aperture photome- 
try with the aperture correction of x 1.33 (§ 2.1.2). Also, we apply 
bootstrap re-sampling (500 times) to obtain the 1-a error-bars for 
the median 24/im fluxes. Note that we include all the Hq emit- 
ters for the stacking, regardless of their individual 24pm detec- 
tion, while we do not use photo-z sources (without Ha detec- 
tion) in order to avoid contaminant galaxies. We derive the total 
infrare d luminosity (Lir) from the stacked 24^m photometry us- 
ing the ICharv & Elbaj ( 1200 ll) SED templates, and then compute 
the SFR using t he com bined Hq and Lir approach suggested by 
iKennicutt et all 120091 ); SFR= 7.9 x 10"^^ [L(HQ)obs + 0.0024 x 
Ltir] (ergs~^). This equation is derived from a tight correlation 
between the combined Hq a nd total IR luminosit ies and the at- 
tenuation corrected SFR (see iKermicutt et al]|2009l) . We note that 
our conclusion does not change even if we calculat e SFR s directly 
from the total IR luminosity using the iKennicut 3 i 19981) relation 
(in this case we tend to derive SFRs which are higher by a fac- 
tor of <2), but the combined Ha and IR approach should be more 
reliable, considering the moderate levels of star-formation activity 
(with ;£10Mq yr~^) of our galaxy sample. 

2.2 Hq emitter samples from our previous studies 

2.2.1 Cluster galaxy sample from MAHALO-Subaru 

In addition to the CI 0939 cluster {z — 0.4), we also use our sim- 
ilar, Hq selected galaxies in a z = 0.8 cluster (RXJ17 16-1-6708) 
and in a z = 2 .2 proto-cluster (PK SI 138— 262), from our previ- 
ous studies (see lKovama et alj20ICl ; 2013 for detailed descriptions 
of these data and the selection of Hq emitters). The Ha emitter 
samples are constructed with narrow-band surveys using NBI19 
(Ac = LlO/im for z = 0.8) and NB2071 (Ac = 2.07Atm for 
z = 2.2) on MOIRCS/Subaru ( ilchikawa et alj|2006l ; ISuzuld et alj 
l2008h . as a part of the MApping H-Alpha and Lines of Oxygen with 
Subaru project (MAHALO-Subaru; see overview bv kodama et al.l 
^013). We use these cluster galaxy samples to study the redshift 
evolution of the SFR-i\f, relation in cluster environments (see 
§ 3.4). We note that our advantage is a perfect matching of their 
redshifts with our similar Ha-selected control field galaxy samples 
from HiZELS (see below). Using only one cluster at each redshift 
might be too simplistic (e.g. our proto-cluster at 2 = 2.2 may not 
necessarily be the progenitor of our z — 0.4 and 0.8 clusters), but 
in this pioneering study, we assume that these Hq emitter samples 
in these three cluster fields represent typical star forming galaxies 
in high-density environments at each redshift. 



2.2.2 Field galaxy sample from HiZELS 

In order to test the environmental dependence of the SFR versus 
Al, relation, we also need a control field galaxy sample. The High- 



© 0000 RAS, MNRAS 000, 000-000 



4 Y. Koyama et al. 



1 r 



-2 



—\ — ' — <- 

z=0.4 



21 



22 23 
z' [mag] 



24 



n 1 1 1 — 

2 2 E- z=0.8 
\. 1 r 





-1 r 



o 



— 1 1 r- 



— 1 1 r- 



-I I I l_ 



_l I l_ 



_l I l_ 



21 



22 
J [mag] 



23 



-| — 1 — 1 — I- 

2^ z=2.2 
1 



-1 

-a E- 



— 1 1 1 r- 



—I I I i_ 



21 



22 23 
K [mag] 



24 



Figure 1. A compaiison between the stellar mass estimates for HiZELS Ha 
emitters derived from the rest-frame S-band magnitudes with one-colour 
correction (see text) and those from full SED fitting derived by Sobral et al. 
(in prep). The line-connected black points show the running median (and its 
associated scatter). It is clear that these two measurements agree with each 
other reasonably well over a wide luminosity range. 



Z Em ission Line Survey (HiZELS: Isest et alJlioTol : ISobral et"ZI 
l2013h is ideally suited for this purpose. This is not only be- 
cause HiZELS is currently the largest narrow-band Ha survey ever 
published, but also because three of their four targeted redshifts 
(z — 0.4/0.8/2.2) are perfectly matched with our cluster sam- 
ples, allowing a direct cluster-field comparison based on the purely 
Ha-selected galaxies. The HiZELS Ha emitter samples are se- 
lected from the UPS and COSMOS fi elds (total ~2 deg^; see also 
iGeach et al.ll2008l : ISobral et al.ll2009l) by wide-field narrow-band 
imaging observations with NB921 on Subaru (2 = 0.4), NBj 
on UKIRT (z = 0.8), and the NBy /H2 filters on UKIRT/VLT 
(z — 2.2; see also iGeach et alj|20l|)^ We select Ha emitters in 
exactly the same way as described in Sobral et^ (l2013h . We note 
that the selection and photometry of the HiZELS sources has been 
made w ith 3" and 2" aper ture for z — 0.4 and z = 0.8/2.2 galax- 
ies (see lSobral et al.ll2013h . respectively, while physical quantities 
of our cluster galaxy samples are measured with total magnitudes 
(MAG_AUTO from SEXTRACTOR). Therefore, we apply an aperture 
correction for the HiZELS sample based on the median difference 
between aperture magnitudes and total magnitudes for each red- 
shift slice, although these corrections are negligibly small in our 
discussion (0.2 dex at maximum). 



2.3 Stellar mass and Ha-based star formation rate 

In this subsection, we derive stellar masses (A/*) and star for- 
mation rates (SFR) of the Ha emitters. The stellar masses of 
galaxies are ideally derived using a SED fitting approach includ- 
ing rest-frame near-infrared (NIR) photometry. However, the rest- 
frame NIR photometry is not available for our cluster galaxy sam- 
ples. We therefore decide to estimate the stellar masses of galax- 
ies based on their rest-frame i?-band magnitudes. The conversions 
from observed magnitudes to M, are determined using the model 



galaxies de veloped by Kodama, Bell, & Bower (see also 

iKodama & A ri motg . 1997 : ^odama et aL.1998i ). and they are ex- 
pressed by the following equations: 

log(M./lo"M0),^o.4 = -0.4(2' -20.07) + A log M0.4, (1) 

log(M./lO"M0),=o.8 = -0.4(J-21.14) + AlogA./o.8, (2) 

log(M./lo"M0),^2.2 = -0.4(i4:(,) - 22.24) + AlogAf2.2. 

(3) 

The final term in each equation (AM) accounts for the colour de- 
pendence of the mass-to-light ratio (M/L) predicted by the same 
galaxy model, and they correspond to: 

A log MoA = 0.054 - 3.81 x exp[-1.28 x (B - z')], (4) 
A log A/0.8 = 0.085 - 2.48 x exp[-1.29 x (i? - J)], (5) 

Alog A/2.2 = 0.030 - 1.50 X exp[-l.ll x {z - K^s))]- (6) 

The stellar masses derived with this "one-colour method" could be 
less accurate compared to those from full SED fitting, but we find 
that there is no systematic difference between the A/, from this 
method and those from the SED-fitting approach using HiZELS 
Ha emitter samples for which full SED information are available 
(see Fig. 1). The two measurements are consistent with each other 
over a wide luminosity range (with a ~0.3 dex scatter), verifying 
that our one-colour method works reasonably well. In this paper, in 
order to make a fair comparison between cluster and field samples 
at different redshifts, we use the above one-colour method for all 
Ha emitters (including MAHALO and HiZELS samples). 

We then calculate SFRs of Ha emitters based on their 
Ha luminosities. We first correct for the contribution of [Nil] 
lines to the total NB fluxes, using an empirical correlation be- 
tween EWre B t(Ha+ [Nll]) and the [Nll]/Ha ratio as described in 
ISobral et al.l ( l2013h . This relation has a large intrinsic scatter, but 
this method is believed to be more realistic approach compared 
with the conventional constant [Nil] correction. We also apply 
a dust extinction correction to the Ha flux of individua l galax - 
ies based on their stellar mass, as shown by iGam & Best! (|2010|) . 
This extinction correction could also be uncertain, given the large 
intrinsic scatter of the Ahq-A/, relation for local galaxies (see 
iGam & Bes3l2010 1. However, it is one of the most reliable esti- 
mators of dust extinction, a nd also the relatio n is reported to be 
unchanged out to z ~ 1.5 dSobral et al.ll2012l) . We therefore ap- 
ply the same extinction correction to all Ha emitters . Final ly, we 
compute the SFRs of galaxies using the iKennicutil 1 1998 1) rela- 
tion: SFRhq(A/0 yr~^)= 7.9 x lO'^^LHcCergs"^). We note that 
the [Nil] line and the dust extinction corrections could be major 
sources of uncertainty, but this is currently inevitable because it is 
impossible to measure [Nil] contribution or dust extinction for in- 
dividual galaxies. We stress that we adopt the same procedure for 
all our sample (cluster and field galaxies), so that we consider the 
relative cluster-field comparison should be robust (see related dis- 
cussions in § 4.5). 

The stellar masses and SFRs derived above may not be very 
accurate for galaxies hosting active galactic nuclei (AGN), as their 
continuum light and Ha line flu xes are contaminat ed b y AGNs. For 
the fie ld Ha emitters (HiZELS). Gam et al.l ( l201oh and lsobral etaTI 
( l2013h carried out a detailed study on the contribution of AGNs into 
their Ha emitters sample. By applying various techniques for AGN 
selection (e.g. X-ray, radio, or emission-line ratios), they find that 
the AGN fraction amongst their Ha emitters is as large as 10% out 
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Figure 2. (Left): The 2-D map of galaxies around the cluster CI 0939 {z = 0.41). Galaxies within both Suprime-Cam and MIPS FoVs are shown. We plot 
24/im-detected Ha emitters, 24/xm-detected photo-2 members (0.30<Zpijot^0-45), all Ha emitters, and all photo-2 members (the meanings of the symbols 
are shown in the plot). We also show the locations of the suiTounding groups following the definitions in Koyama et al. (201 1). Contours are drawn based on 
the surface number density of all member galaxies (same as fig. 5 of Koyama et al. 2011). (Right): The colour-magnitude diagram for each environment as 
defined in the left-hand panel. The meanings of the symbols are the same as the left-hand panel. The (blue) hatched and (orange) cross-hatched histograms 
show the colour distribution of Ha emitters and 24/im sources, respectively. 



to z lor~15%at2 > 1. For cluster samples, we confirmed a 
few (up to five) AGNs in each cluster using X-ray imaging data or 
spectroscopic information (see Koyama et al. 2008; 2011; 2013). 
Considering the size of the total galaxy sample we use in this paper 
(~100 galaxies per cluster), it is roughly estimated that the AGN 
fraction in our cluster galaxies is ~5-10%. Unfortunately, it is not 
possible to fully quantify the contribution of AGNs in our cluster 
Ha emitters in the same way as for HiZELS samples (due to the 
lack of multi-wavelength data). In order to do a fair comparison 
between cluster and field samples, we do not exclude AGNs from 
our cluster or field Ha emitters in the following discussions, but we 
note again that the AGN contribution is always small (~10%) and 
should not strongly bias the results. 



3 RESULTS AND DISCUSSION 

3.1 Panoramic Ha and MIR view of the CI 0939 cluster 

In Fig. 2, we show the spatial distribution of the MIR-detected Ha 
emitters and the MIR-detected photo-z member galaxies. We also 
show all the Ha emitters and photo-2: member galaxies, and the 
locations of the West Clump and the fo ur surr ounding groups fol- 
lowing the definitions in Kov ama et alj feoill) . Note that we only 
show the galaxies located within the overlapped regions between 
our Subaru and MIPS FoVs. In the right-hand panel of Fig. 2, we 
show the colour-magnitude diagram for cluster (1 Mpc from the 
cluster centre), group (including the west clump and four groups), 
and field environments. These colour-magnitude diagrams show an 
overall trend that the MlPS-detected sources tend to be luminous 
(z' < 21 mag) and to have red optical colours (B—I > 2 mag). 

We compare the B—I colour distribution of the Ha emitters in 
each environment (see histograms in Fig. 2). The lack of the red Ha 
emitters (B—l>2 mag) is clearly visible in the cluster environment. 



while a large number of Ha emitters tend to have red colours in 
the group environment. The Kolmogorov-Smimov (KS) test shows 
that the probability that the colour distributions of Ha emitters in 
clusters and groups environment are from the same parent popula- 
tion is 1%. The trend becomes less significant if we compare the 
group galaxies with field galaxies, but the KS test still suggests that 
the group and field galaxies are unlikely to be drawn from the same 
parent population (2.5%). This probability goes down to <2% if 
we use "all" Ha emitters sample within the Suprime-Cam FoV 
(including the galaxies located outside the MIPS data coverage). 
Furthermore, we calculate the fraction of the red galaxies with any 
indication of star formation activity (i.e. Ha or MIR detection): 
/=A^rodSF/A^anrod- We find fractions of 4±1%, 17±3%, 11±1% 
for cluster, group, and field environments, respectively. We there- 
fore conclude that the red star-forming galaxies are most frequently 
seen in the gro up-scale env ironment (supporting our previous find- 
ina in Kovama etaLll201lh ; in other words, a non-negligible frac- 
tion of optically red galaxies in the group environment are still ac- 
tively forming stars. 



3.2 Stacked MIR properties as a function of galaxy colours 

The colour-magnitude diagrams in Fig. 2 show that a fraction of 
the red Ha emitters are individually detected at 24/^m, suggesting 
they are dusty sources rather than passive galaxies. However, the 
limited depth of the 24/im data prevents us from assessing the gen- 
eral MIR properties of the faint Ha galaxies. We therefore apply a 
24/im stacking analysis to study the properties of the Ha emitters 
more generally. To do this, we divide the full Ha emitter sample in 
the CI 0939 field into three colour bins (at B-I=\.5 and 2.0), and 
perform stacking analysis as described in § 2.1.3. The results are 
shown in Fig. 3a. It is clear that the redder Ha sources have higher 
SFRs than bluer sources. We also show in Fig. 3a the SFRs derived 
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Figure 3. The stacked properties of the Ha emitters as a function of B—/ colours (rest-frame U—V). The vertical error-bars show the Icr uncertainty derived 
from the bootstrap re-sampling approach during the MIPS stacking (see text), while the horizontal eiTor-bars show the 25-75% colour distribution of each 
subsample. (a) The SFRs derived from 24^m stacking analysis (black symbols). The blue symbols indicate the median values of SFRjjq (without dust 
extinction correction) for each subsample. (b) The Ajja value for each subsample, calculated from the ratio of SFRjjs and SFRhq. (c) The median stellar 
mass for each colour subsample (black), as well as those for individual sources (grey dots), (d) The SSFR for each subsample derived as SFR nonnalized by 
stellar mass. These four plots demonstrate that the red Ha emitters are massive star-forming galaxies with higher SFR and higher dust extinction compared 
with normal blue Hq emitters. 



from the median Hq flux for each subsample, without dust extinc- 
tion correction. The same trend is still visible, but it is weaker than 
that obtained from MIPS stacking analysis. This result suggests that 
the redder Ha galaxies are dustier than bluer galaxies. In Fig. 3b, 
we show the extinction at Ha derived from a ratio of the total SFR 
(from MIPS stacking) and the dust-free Hq SFR as a function of 
galaxy colours. This plot clearly shows that the red Ha emitters 
have much higher extinction (with Ahq^LS mag) compared with 
the blue emitters (ylna^l mag), suggesting that the red Ha emit- 
ters are dusty, star-forming galaxies. 

We should note that the clear trend that the redder Hq emit- 
ters have higher SFRs could be produced by a stellar mass differ- 
ence between red and blue Ha emitters. In Fig. 3c, we show the 
estimated M, for the individual Hq emitters (see § 2.3), as well 
as the median stellar mass for each colour subsample. The red Ha 
emitters tend to have much higher stellar masses (by a factor of 
~10) than the blue Hq emitters, which is more significant than 
the difference in SFR (i.e. Fig. 3a). Consequently, it turns out that 
the red Hq emitters tend to have lower specific star formation rate 
(SSFR) than blue Hq galaxies (see Fig. 3d). However, the SSFR 
for the red Hq emitters are still at the ~10~^" yr~^ level, which 
is about one order of magnitude higher than the usual definitions 
of "passive" galaxies (10~^^ yr~^), suggesting they are part of the 
"star-forming" population. 

As a fuither check, we show in Fig. 4 (left) the results of our 
stacking results for the red and blue Hq emitters on the SFR-Af, 
diagram (we divide the red and blue Ha emitter sample into two 
equal-sized stellar mass bins). Our data points show an excellent 
agreement with the best -fitted SFR-M<. rel ation for the z — OA Ha 
emitters from HiZELS ( Sobral et al. I l2013 h. further supporting our 
conclusion that the red Hq emitters are dusty star-forming galax- 
ies, rather than dust-free passive galaxies. We recall that when de- 
riving SFR- A/* relation for the HiZELS sample, we adopt an Af,- 
dependent dust extinction correction to HiZELS data (see § 2.3). 
The excellent agreement (over the wide Af, range) between the 
SFRs derived from the MIPS stacking analysis and those from the 
independent, extinction corrected Hq fluxes also supports the va- 



lidity of our procedure for the dust extinction correction applied in 
this study. 

3.3 Environmental dependence 

We have shown that our Ha emitters tend to be located on the gen- 
eral SFR-Af, sequence (Fig 4-left). An interesting question here 
is: does environment influence the SFR- A/* relation? To answer 
this question, we attempt a similar analysis to § 3.2 but dividing the 
whole 2 = 0.4 Hq emitter sample into four environment bins based 
on the local galaxy density (at log Eio =1.6, 1.9, 2.2). The density 
is calculated using all cluster member galaxies (photo- z selected 
and Hq selected) with the nearest-neighbour approach, calculated 
within a radius to the lOth-nearest neighbour from each source. We 
further divide each environment subsample into two equal-sized 
stellar mass bins, and perform the 24fim stacking analysis. The re- 
sults are shown in Fig. 4 (right). The stacking analysis becomes 
challenging, particularly for the lower-mass sources due to the lim- 
ited sample size, but broadly speaking, all our data points are likely 
to be located on the same SFR-Af, sequence. This suggests that the 
SFR-Af, (or SSFR-Af,) relation for star-forming galaxies does not 
strongly correlate with the environment at z = 0.4. 

Using the same environmental subsamples, we also test the 
SFR-density relation for star-forming galaxies at z — 0.4. In Fig. 5 
(left), we show the SFRs from 24//m stacking analysis as a func- 
tion of galaxy density. It is interesting to note that the SFR increases 
with increasing galaxy density (by a factor of ~3-4), and the trend 
is confirmed for both red and blue Ha emitter samples. In contrast, 
this trend becomes much less significant when we normalise the 
SFRs by Af» to compare the specific SFR (see Fig. 5 right). The 
fact that the SSFR shows a weaker environmental trend than SFR 
(Fig. 5) implies that the SFR excess detected in the high-density 
environment would largely be explained by the M, difference be- 
tween the different environments. We note, however, that the Af, 
distribution amongst Ha emitters does not seem to be strongly de- 
pendend on environment (as shown by histograms in Fig. 4b). We 
can still find a weak trend that the galaxies in high-density environ- 
ment tend to be more massive, but the difference is at the ~0.2 dex 
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Figure 4. (Left): The median SFRs from 24fim stacking analysis for the red and blue Ha emitters as a function of stellar mass. The red and blue samples 
are further divided into two equal-sized stellar mass bins. The dotted line and the shaded region show the best-fitted SFR-Af, relation and its scatter for the 
HiZELS z = 0.4 Ho emitter sample from Sobral et al. (2013). All our data points show an excellent agreement with the HiZELS relation. (Right): The same 
plot as the left-hand panel, but for our environmental subsamples (each sample is further divided into two equal-sized stellar-mass bins). In both panels, we 
show the A/* distribution (histograms) and the median stellar masses (aiTows) for each sample. In each plot, the vertical error-bars are from the bootstrap 
re-sampling in the MIPS stacking analysis, while the horizontal error-bars show the 25-75% A/, distribution for each subsample. 



level at maximum. The weakness or lack of the environmental 
dependence of the stellar mass distribution a mongst star-forming 
galaxies is consiste nt with some recent studies dOreene et al.l2013 : 
iGiodini et al.ll2012h . but this small Af» difference may not be able 
to fully account for the significant SFR increase towards high- 
density environment. We therefore speculate that the SFR excess in 
the high-density regions could be explained by a "mixed effect" of 
both slightly higher stellar masses and a small SSFR excess (both at 
~0.2 dex level) in the high-density environment. Indeed, the small 
SSFR excess is visible as the remaining positive slope in the SSFR 
versus Eio plot (see Fig. 5 right), which is also equivalent to the 
small positive offset of the stacked data points for high-density 
environments in the SFR-A/, plot (Fig. 4 right). However, a firm 
conclusion will require future studies with larger galaxy samples, 
and so our reasonable conclusion at this stage is that the overall 
SFR-M* relation is not strongly correlated with the environment 
at 2: = 0.4. 



It should be noted again that the star-forming activity of galax- 
ies presented here only focuses on star-forming galaxies, and that 
we do not include passive gala xies in our analysis. Indeed, we 
showed in iKovama et all llOll i) that the Ha emitter fraction is a 
strong function of environment, showing a significant decline to- 
wards the cluster core. The important message from our current 
analysis is that the SFRs of the uniformly Ha-selected galaxies 
depend on enviroimient in the sense that galaxies in high-density 
environment have higher SFRs. This would suggest that the star- 
forming galaxies "surviving" in the very high-density environment 
have experienced an accelerated stellar mass growth in their past 
(e.g. through starburst activity), while it is also possible that they 
have been forming stars at constant rates over a longer period of 
time than general field galaxies (i.e. nature effects). Unfortunately, 
we cannot discuss the past of their star formation history with our 
photometric data alone. In a future paper, we will investigate their 
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Figure 5. The SFR (left) and SSFR (right) for our Ha emitters derived 
from 24/^m stacking as a function of the local galaxy density (Sio). The 
density is calculated with all member galaxies (photo-z members and Ha 
emitters), while we include only Ha emitters in the stacking analysis. The 
SFR increases with increasing density regardless of their colours, while the 
trend becomes much flatter if we compare SSFR. The vertical error-bars 
show the 1-(T distribution from the bootstrap resampling and the horizontal 
error-bars show the 25-75% distribution, same as previous figures. 



spectroscopic properties through an intensive spectroscopic follow 
up programme on this CI 0939 field (Nakata et al. in preparation). 



3.4 The evolving SFR-M, relation since z~2 

We have examined the star-forming activity of galaxies and its de- 
pendence on stellar mass and environment at 2 = 0.4. Here, we dis- 
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Figure 6. The SFR-A/» relation for the Ha-selected galaxies at z = 0.4, 0.8, 2.2. In each panel, we plot the Hcj emitters in cluster environment located 
within Rc <1 Mpc and at 1< Rc <2 Mpc from the cluster centre. We also plot the Ha emitters in general field environment selected from HiZELS at 
the same redshifts. Note that all the Ha emitters plotted here are selected as those having EW(HQ+[NlI])rcst>30A. The dotted lines show the SFR and Af, 
cut applied at each redshift. The dashed line shows the best-fitted SFR-Af» relation for the HiZELS sample, while the solid line shows the relation for the 
cluster sample assuming the same slope as the HiZELS relation. The dot-dashed fine is the local (z = 0) relation derived from the equation provided by 
Whitaker et al. (2012). It is clear that the SFR-A/* relation evolves significantly since z ~ 2 in both cluster and field environment, while at fixed redshifts, 
the environmental dependence of the SFR-Af, relation seems to be very small. 



cuss the environmental dependence of the evolution of star-forming 
galaxies in a broader context, particularly focusing on the evolu- 
tion and environmental dependence of the SFR versus A/, relation 
across cosmic time. By compiling all of the Ha emitter samples to- 
gether (including our MAHALO and HiZELS samples; see § 2.2), 
we examine the SFR versus A/» sequence in cluster and field envi- 
ronments al z = 0.4, 0.8, and 2.2. We note that, while all the sam- 
ples are selected based on the Ha line, the EW cut applied in each 
survey is slightly different. Therefore we decide to apply the same 
(rest-frame) EW cut to all the sample, EWrcst(Ha-l-[Nll])=30A, 
down to which all our Ha data are complete. We also note that it is 
not possible to make a strictly fair comparison between clusters at 
different redshifts; for example, a density-based definition of envi- 
ronment requires accurate membership determination, while using 
the virial radius could be misleading because our high-z clusters are 
not virialised yet. Therefore, we instead use galaxies within 2 Mpc 
(in physical scale) from each cluster centre as "cluster" galaxies in 
the following discussion. 

In Fig. 6, we show all the Ha-selected galaxies in clusters and 
field environments at each redshift. We use different symbols for 
the cluster Ha emitters located within i?c <1 Mpc, l<_Rc<2 Mpc, 
and the field Ha emitters from HiZELS (see labels in the plot). For 
z = 0.4, we also show the Ha emitters located at Rc >2 Mpc from 
the CI 0939 cluster. This plot clearly shows that the SFR- A/, rela- 
tion evolves with redshift, while the relation is always independent 
of environment out to z ~ 2, q ualitatively consistent with the situ- 
ation in the local Universe (e.g. IPeng et al]|2010l ; IWiiesinghe et al.l 
I2OI2I) . We also show the best-fitted SFR-Af, relation in each panel. 
It can be seen that the offset between the relation for cluster and 
field galaxies is always small. The slope of the SFR-A/* relation 



tends to be steeper for the lower-redshift samples, which is also 
qualitatively consistent with previous works (but we note that the 
best-fit relation drawn on the plot is uncertain because our star- 
forming galaxy samples are not completely stellar-mass limited). 

We here examine the galaxy distribution on the SFR-Af, 
plane more in detail. In Fig. 7, we show the distributions of A/,, 
SFR, and the offset from the main sequence (best-fitted SFR-A/, 
for field galaxies). The shaded histograms are for field Ha emitters, 
and the hatched histograms are for cluster Ha emitters. Broadly 
speaking, it seems that there is no significant environmental dif- 
ference at any of the three epochs. For z = 0.4/0.8, the KS test 
suggests that it is unlikely that the cluster and field galaxies are 
from a different parent population, while for z = 2.2 we find 
a possible trend that cluster galaxies have a small excess in all 
three properties (KS test shows <1% probability that the clus- 
ter and field samples are from the same parent population). We 
note that this trend is qualitatively consistent with some earlier 
studies showing a higher stellar mass in star-forming galaxies in 
z > 2 proto-clij sters (e.g. Steidel et al. 2005; Kurketal. 200,1; 
iHatch et alj[20Tll ; iMatsuda et alJi201 ll ; iKovama et alil2013i ). Our 
current analysis also supports the idea that the star-forming galax- 
ies in proto-cluster environment tend to be more massive than the 
general field galaxies, and this Af, excess would also account for 
the SFR excess in the proto-cluster environment; this may repre- 
sent an accelerated galaxy growth in the early phase of the cluster 
assembly history. 

Finally, we quantify the evolution of the star-forming activ- 
ity of star-forming galaxies in the cluster and field environment. In 
Fig. 8, we plot the redshift evolution of the SSFR of Ha emitters 
in clusters (red squares) and in field environment (black squares) at 
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the stellar mass of log{M*/MQ) = 10. It is found that the SSFR of 
Ha-selected galaxies evolves significantly, going approximately as 
(1 + z)^, since z ~ 2 in both clusters and the field. This decline 
of SSFR is in good agree ment with many studies of the cosmi c 
star formation history (e.g. lYoshida et al .l2006tlKarim etalj201l[) . 
Therefore, an important indication from this study is that the evo- 
lution of star-forming galaxies in cluster environments just follows 
the same evolutionary track as that of general field galaxies. 



3.5 Comparison with other studies 

Studying the environmental dependence of galaxy star formation 
activity in the distant Universe is obviously an important step to- 
wards understanding the physical processes which drive the envi- 
ronmental effects. Since discus sion on the "revers al" of the SFR- 
density relation was invoked bv lElbaz et al. I( l2007h . there has been 
much debate about the role of environment in the distant Universe. 
We note that our results ar e qualitatively consi stent with some re- 
cent studies. For example, ISobral et al.l (1201 ih used Ha emitters 
sample sX z = 0.8 selected from HiZELS to show higher median 
SFRs (by a factor of ^--^2-3) in high-density environment compared 
to low-density environment. They also showed that the stellar mass 
of Ha emitters is weakly correlated with the environment (with 
~0.3 dex increas e in their h i ghest- density bins), which also agrees 
with our finding. iTran et al] ( l2009t) analysed MIR data of a "super- 
group" environment at z = 0.37 to show that the characteristic IR 
luminosity (L*) in the group environment is higher than that in the 
field, based on the analysis of the IR luminosity function (see also 
IChung et alj2oT(]|) . Related to this, some studies of distant clusters 
show a peak of star formation activity at a certain galaxy density 
which corresponds to group or cluster outskirts envirormien t (e.g. 
IPoggianti et ai]|2008l ; [Kovama et alj2oicl : lGeach et aLll201 lb . 

On the other hand, a growing number of studies recently have 
reported a weakness or absence of any relation between SSFR and 
environmental density amongst star-forming galaxies at least out to 
z 1 , or possibly to z ~ 2. In the local Universe. iBalogh et al.l 
( I2OO4 ') showed that the EW(Hq) distribution (equivalent to SSFR 
distribution) amongst star-forming galaxies is independent of en- 
vironment. More recent studies also indicated that the SFR-M, 
relation for local star-forming galaxies does not correlate w ith the 
environment (e.g. lPeng et alj|2010l : lwiiesinghe et al.ll2012h . Simi- 
lar suggestions have also been made fo r distan t star-forming galax- 
ies as well. For example. IMcGee et al. J) studied a large sam- 
ple of z = 0.4 group galaxies to show that the average SSFRs of 
star-fo rming galaxies are t he same in groups as in field environ- 
ments. lMuzzinetal]j2012 ) also showed that SSFR of star-forming 
galaxies is independent of environment at fixed stellar mass from 
their detailed spectrosco pic survey of 2 ~ 1 cluster galaxies (see 
also iGreene et al.ll2012h . Furthermore, our recent studies of dis- 
tant (proto-)clusters also find a hint that the SSFR of star-forming 
galaxies is indepe ndent of envirormient at fixed stellar m ass out 
to z ^ 2.5 (e.g. iTadaki et alj|2012l : [Havashi et alj|201 ll : 2012; 
iKovama et alj2013[K 

However, it should be noted that the independence of the SFR- 
Mt relation for star-forming galaxies with environment i n the dis- 
tant Universe is still controversial (Patel et al. 2011; Vulc ani et al.l 
l20ld ;l Li et al.ll201ll) . Indeed, this kind of analysis could be highly 
sensitive to the sample selection, the measuremen t of star formatio n 
rates, and the definitions of environment (see e.g. IPatel et al.l2009h . 
Our data are purely Ha selected (for both cluster and field galax- 
ies), and in this sense our cluster-field comparison would be robust. 
One possible bias is that our data are complete only for relatively 
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Figure 7. The distribution of M* (left), SFR (middle), and offsets from 
the main sequence of field galaxies (right) at each redshift. The shaded his- 
tograms show the results for HiZELS sample, while the hatched histograms 
show the results for cluster {Rc < 2 Mpc) galaxies. The vertical dotted 
lines in the left and middle panels show the A/» or SFR cut we applied for 
each redshift sample, while the dotted lines in the right-hand panel show the 
location of the zero-offset. The actual difference between cluster and field 
galaxies is always small (^^0. 1-0.2 dex at maximum), but we note that a 
statistical test suggests that the two distributions may be different for our 
z = 2.2 sample in the sense that the cluster galaxies have higher Af» and 
higher SFR (see text). 



strong emitters; we recall that our definition of star-forming galax- 
ies is EW>30A, so that we cannot discuss faint, low-EW sources. 
It is likely that such low EW sources do exist in both environments 
(and perhaps they may be more numerous i n cluster environment). 
However, as reported in ISobraletal.l(l20T 1]), such low-EW sources 
tend to be dominated by massive galaxies with relatively low SFR 
(i.e. largely "switched-ofT' population), so that we expect that it is 
unlikely that such low-EW sources would have a significant impact 
on our discussions on "star-forming" galaxies. 



3.6 Interpretation and Caveats 

Our main finding of this study is that the SFR versus A/, relation 
for the Ha-selected galaxies does not strongly depend on the en- 
vironment at any time in the h istory of the Univ erse since 2 ~ 2. 
This is a similar suggestion bv |Pengetal1 ( l20IOl) . who used COS- 
MOS data to study the redshift evolution of the SSFR of (blue) 
star-forming galaxies, and find no environmental difference since 
z ~ 1. The independence of the SFR-A/» relation with envi- 
ronment could be ex plained if the envir onmental quenching is a 
rapid process (see e.g. lMuzzin et al .l2012h . That is, the environment 
instantly shuts down the star-formation activity of galaxies once 
the environmental effects are switched on, so that declined star- 
formation is not observed (because our galaxy samples are selected 
with Ha). Therefore a naive interpretation of our result would be 
that the major environment quenching mechanisms are always fast- 
acting in the history of the Universe since z ~ 2. 

An important, but unexplored issue is the contribution of 
AGNs. While most of our Ha-selected galaxies are likely to be 
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Figure 8. The redshift evolution of the SSFR at M, = IO^^A/q de- 
lived from the best-fitted SFR-Af» relation for cluster (red squares) and 
field (black circles) galaxies. The error-bars represent the standard deviation 
around their best-fit relation (see Fig. 6). The dotted lines are the evolution- 
ary tracks following oc (1 + z)^, oc (1 + 2)'^, and oc (1 + 2)^, to guide 
the eye. The local data point is derived by adopting 2 = in the equation 
ofWhitaker et al. (2012). 



powered by star formation (see § 2.3), there still remains a possi- 
bility tiiat the AGN contribution coul d be dependent upon redshift, 
mass, and environment. For example. IPopesso et al ] (|20T ij) carried 
out a detailed FIR study of the star forming activity of galaxies at 
z ^ 1 using Herschel data. They find that, while overall the SSFR- 
A/* relation does not depend on environment, the reversal of the 
SFR-density relation could be produced by very massive galaxy 
population. They also noted that the inclusion of AGNs into the 
analysis could also lead to an apparent reversal of the SFR-density 
relation. Therefore, more detailed studies of individual galaxies 
(including spectroscopy) are clearly needed to unveil the role of 
AGNs, as a future step of this study. 

Another caveat on our result concerns the prediction of dust 
extinction correction. We applied the empirical correction based on 
the AYia-Mf correlation establishe d for local galaxies (see § 2.3), 
which has a large intrinsic scatter (iGam & Be st"2010l). The rela- 
tion is reported to be unchanged out to z ~ 1.5 ( Garn et al. 20 13; 
ISobraletal]|2012l) . and as shown in ICam et al.l ( |2010|) . the A^a 
values (from the ratio of IR and Ha luminosities) show no sig- 
nificant environmental variations. Indeed, as we showed in § 3.3 
(Fig. 4), our results are unchanged even if we derive SFRs from a 
MIR stacking analysis, suggesting that this effect is not large. Us- 
ing our CI 0939 sample presented in this paper, we further check 
the extinction values {Ana from MIPS stacking analysis) for clus- 
ter {Ra <2 Mpc) and field (i?c >2 Mpc) galaxies. We find that 
the field galaxies have slightly (~0.2 mag) higher extinction, while 
we find an opposite trend if we use the density-coded definitions of 
environment (galaxies in high-density environment tend to be more 
extinguished). Therefore, it is not possible to say anything about the 
environmental dependence of the dust properties of galaxies from 
our data alone (due to limited sample size), but in any case, it is 
unlikely that the dust properties in star-forming galaxies are signif- 
icantly different between different environments. 



4 SUMMARY AND CONCLUSIONS 

In this paper, we study the evolution and environmental dependence 
of the SFR-A/, correlation for star-forming galaxies since 2 ~ 2. 
We first present the MIR properties of the Ha-selected star-forming 
galaxies in a rich cluster at 2 = 0.4 (CI 0939), and then we com- 
pare the z = 0.4 galaxies with our similar, Ha-selected galaxies at 
different redshifts and in different environments. Our findings are 
summarized as follows: 

(1) The red Ha emitters, which are reported to be most fre- 
quently seen in the group-scale environment at 2 = 0.4, are dusty 
red galaxies rather than (semi-)passive galaxies. Using a wide-field 
SpitzerRs/SPS 24/im dataset, we find that a large number of massive 
red Ha sources are individually detected at 24/im, suggesting they 
are luminous and dusty. Also, with a stacking analysis, we confirm 
a more general trend that the red Ha sources tend to have higher 
SFRs with stronger dust extinction compared with normal blue Ha 
emitters. 

(2) We also find that the median SFR of Ha emitters (derived 
from the MIR stacking analysis) increases with increasing galaxy 
number density at 2 = 0.4. This result is confirmed for both red 
and blue Ha emitters, while the trend becomes much weaker if we 
compare their SSFR. We note that there still remains a positive cor- 
relation between SSFR and galaxy number density, and therefore 
we speculate that the SFR excess in the high-density environment 
can be caused by a mixed effect of both slightly higher A/* and a 
small SSFR excess (both at ^--^0.2 dex level) in high-density envi- 
ronment, implying that the Af, difference may not be able to fully 
account for the SFR increase in the high-density regions. 

(3) Using our large Ha-selected galaxy samples in distant 
cluster environments (from MAHALO-Subaru) and in general field 
environments (from HiZELS) at z — 0.4, 0.8, 2.2, we examine 
the environmental dependence of the SFR- A/* relation across cos- 
mic time. We find that the SFR-A/* relation evolves with cosmic 
time, but that there is no detectable environmental variation in the 
SFR-Af, relation at any of these redshifts. This result suggests 
that, at fixed redshift, the overall star-forming activity amongst star- 
forming galaxies is not strongly correlated with environment once 
their stellar mass is fixed. On the other hand, for the 2 — 2.2 sam- 
ples, we find a possible trend that the (proto-)cluster galaxies tend 
to have higher M, (hence higher SFR), suggesting an accelerated 
galaxy mass growth in the early phase of cluster assembly. 

(4) Based on the SFR-A^* relation we derived for cluster and 
field galaxies at 2 = 0.4, 0.8, 2.2, we also examine the evolution 
of the SSFR for star-forming galaxies (at the fixed mass of A/* = 
lO^^Afo). We find that the SSFR evolves significantly, as (1 + 
2)^, in both cluster and field environments. This means that the 
star-forming galaxy evolution in cluster environments just follows 
the same evolutionary track as that of field galaxies, which is most 
simply inteipreted as implying that the primary physical driver of 
the environmental quenching is always a fast-acting process at any 
time in the history of the Universe since 2 ~ 2. 
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